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ABSTRACT 

We have investigated chemical evolution in the young universe by analysing the detailed chemical 
enrichment pattern of a metal-rich galaxy at high redshift. The recent detection of over 20 elements 
in the gas-phase of a damped Lyman-a absorber (DLA) at z = 2.626 represents an exciting new 
avenue for exploring early nucleosynthesis. Given a strict upper age of ~ 2.5 Gyr and a gas-phase 
metallicity about one third solar, we have shown the DLA abundance pattern to be consistent with 
the predictions of a chemical evolution model in which the interstellar enrichment is dominated by 
massive stars with a small contribution from Type la supernovae. Discrepancies between the empirical 
data and the models are used to highlight outstanding issues in nucleosynthesis theory, including a 
tendency for Type II supernovae models to overestimate the magnitude of the "odd-even" effect at 
subsolar metallicities. Our results suggest a possible need for supplemental sources of magnesium and 
zinc, beyond that provided by massive stars. 

Subject headings: galaxies: abundances — galaxies: chemical evolution — galaxies: intergalactic 
medium 



1. INTRODUCTION 

Chemical abundances in local stars are fossil evidence 
from which early nucleosynthetic processes may be in- 
ferred. An alternative probe into the earliest epochs of 
the universe is via the direct detection of metals in high 
redshift neutral gas. This requires a suitable background 
source such as a quasar, whose spectrum reveals absorp- 
tion features in the intervening gas clouds. Damped 
Lyman-a Absorbers (DLAs) are one such class of quasar 
absorption system, defined as having an HI column den- 
sity N{m) > 2 X 10^" cm~2 

The recent discovery of an intervening galaxy along 
the sightline to the quasar FJ081240.64-320808 (hereafter 
DLA-B/FJ0812+32) has lead to the most comprehensive 
abundance pattern measured beyond the local universe 
(Prochaska, Howk, & Wolfe 2003, hereafter PHW03). The 
gas-phase abundances of over 20 elements - including B, Cu, 
Ga, and Ge - have been detected in this galaxy, which has an 
unusually high column density and metallicity. Many of these 
elements have never before been measured at high redshift 
and provide a unique opportunity to probe early conditions 
of galaxy formation. 

Identifying the modern-day counterparts to DLAs is still an 
open question. Kinematic evidence is consistent with DLAs 
being large protothick disks (Prochaska & Wolfe 1997) or 
smaller merging protogalactic clumps (Haehnelt et al. 1998), 
while chemical properties have lead others to favour inter- 
vening dwarf or irregular galaxies (e.g. Lanfranchi & Friaga 
2003 and references therein) . Galaxies giving rise to DLAs are 
likely to encompass a wide size and morphological range, from 
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dwarf irregulars to giant ellipticals. Indeed, imaging studies 
of low redshift DLAs expose a variety of galactic types (e.g. 
Le Brun et al. 1997; Rao et al. 2003). Observed DLAs may 
tend to probe different types of galaxies at different evolu- 
tionary stages as a function of redshift. 

The metal-line absorption profiles of the galaxy DLA- 
B/FJ0812-I-32 provide a snapshot of the chemical composi- 
tion in the interstellar gas along only one sightline at one 
point in time. However a sufficiently detailed abundance pat- 
tern opens a window into the past history of star formation 
because chemical elements are produced on different charac- 
teristic timescales depending on their nucleosynthetic origin. 
The abundance of metals in the gas-phase of galaxies and pro- 
togalaxies refiects the interplay between processes including 
star formation, the synthesis of elements by stars, the release 
of metal-rich material via stellar winds and supernovae, and 
the fiow of gas. We have constructed a range of simple chemi- 
cal evolution models that follow these processes as a means of 
analysing this DLA. We discuss instances where the inability 
to correctly predict specific elemental abundances may indi- 
cate gaps in our understanding of nucleosynthesis. 

In Section |21 we describe the observations of the DLA at 
z = 2.626 along the sight-line to FJ081240.6-)-320808, upon 
which this study is based. Section [^describes the ingredients 
of the chemical evolution code used to model this DLA. In 
Section |1] the theoretical predictions are compared with the 
DLA data and we discuss the sensitivity of the results to age, 
metallicity, and model ingredients. 

2. OBSERVATIONS 

DLA-B/FJ0812-I-32 was identified in the discovery spec- 
tra of the FIRST Bright Quasar Survey (White et al. 2000). 
Subsequent R ~ 8000 spectroscopy with the ESI spectrom- 
eter (Sheinis et al. 2002) revealed its very strong metal-line 
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Table 1 

Elemental Abundances in DLA-B/FJ0812+32 ° 
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[X/Hf 




5dc(90% c.l. 


r [x/s]- 
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-0.57 


0.085 


0.1 (0.05) 


-fO.3 
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-2.24 


0.058 


0.0 (0.1) 


> -1.47 


O 




-0.54 


0.101 


0.1 (0.05) 


+0.33 


Mg 




-0.78 


0.053 


0.3 (0.1) 


+0.29 


Al 
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-2.00 


0.054 


> 0.5 


> -0.73 


Si 




-0.91 


0.053 


0.3 (0.1) 


+0.16 
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-1.06 


0.000 


< 0.3 


< +0.01 
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-0.87 
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0.1 (0.05) 
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CI 




-1.55 


0.000 


> 0.0 


> -0.78 


Ti 




-1.87 


0.112 


> 0.7 


> -0.4 


Cr 




-1.61 


0.032 


> 0.7 


> -0.14 


Mn 
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-1.85 


0.000 


0.7 (0.1) 


< -0.38 


Fe 




-1.69 


0.017 


> 0.7 


> -0.22 


Co 
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-1.48 


0.000 


> 0.7 


> -0.01 


Ni 
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> 0.7 


> -0.26 


Cu 
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-1.11 


0.000 


> 0.7 


< +0.36 


Zn 
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0.2 (0.1) 


+0.06 


Ga 
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-1.45 


0.000 


0.7 (0.1) 


< +0.02 


Ge 




-0.92 


0.035 


0.3 (0.1) 


+0.15 


As 


< 0.26 


0.000 
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< +1.03 


Kr 
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-0.44 


0.000 


0.0 (0.1) 


< +0.33 


Sn 
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-0.27 


0.000 


0.0 (0.1) 


< +0.5 


Pb 
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-0.10 


0.000 


0.0 (0.1) 


< +0.67 



"Measurements taken by Prochaska, Howk, Sz Wolfe (2003). 
''Gas-phase abundance on a logarithmic scale relative to solar, where 
A^{HI) ^ 102135 cm"^. 

"^Statistical error on gas-phase abundances. 

^Dust-corrections and uncertainties estimated from depletions patterns 
observed in Galactic gas. 

"^Dust-corrected abundances on a logarithmic scale relative to S. 



transitions (Prochaska et al. 2003). PHW03 published HIRES 
echelle spectra of this 'metal-rich' DLA and reported the de- 
tection of over 20 elements. The resulting abundance pattern 
provides an unparalleled laboratory for the study of elemental 
production in the young universe. 

PHW03 derived gas-phase abundances using standard tech- 
niques and estimated uncertainties in the measurements by 
propagating the statistical errors. With only a few excep- 
tions, the dominant uncertainty in the resulting elemental 
abundances is the effect of differential depletion. Similar to 
the ISM of the Milky Way, refractory elements like Ni, Fe, 
Cr may be depleted from the gas-phase such that their ob- 
served abundances are significantly lower than their intrinsic 
values. This is the principal challenge to studying chemical 
abundances in the DLA and it has led to competing inter- 
pretations of their chemical evolution history (e.g. Lu et al. 
1996; Kulkarni, Fall, & Truran 1997; Prochaska & Wolfe 2002; 
Vladilo 2002). In the following analysis, we adopt the conser- 
vative dust corrections imposed by PHW03 based on empiri- 
cal depletion patterns observed for the Milky Way ISM (e.g. 
Savage & Sembach 1996). Table Q lists both the gas- phase 
abundances and dust-corrected abundance ratios for the ele- 
ments observed in DLA-B/FJ0812+32 by PHW03, along with 
statistical errors and dust-correction factors. When relevant, 
we discuss the implications of dust depletion for our conclu- 
sions. Further detail regarding the observational errors will 
be given in Prochaska et al. (2003, in preparation). 

3. THE MODEL 

We simulated the chemical enrichment history of a DLA 
using the chemical evolution code described in Fenner & Gib- 



son (2003). Under this formalism, primordial gas is allowed 
to collect and form stars that synthesise new elements. An 
exponentially decaying gas infall rate was adopted and there 
were no outfiows of gas driven by galactic winds. In this way, 
the models described below resemble different regions of disk 
galaxies. The evolution of the gas phase abundance pattern 
reflects the cumulative history of the dynamic processes of 
star formation, stellar evolution, and nucleosynthesis. The 
classic set of equations governing these processes (described 
in the seminal work of Tinsley 1980), have been numerically 
solved by defining ai{t) as the mass surface density of species 
i at time t, and assuming that the rate of change of ai{t) is 
given by: 



— cr,(t) = / V(i - r^) Y,{m, Z{t - r™)) -^-^ dm 
d . 

+ -^0-i(t)i„fall 

-x^{t)m, (1) 

where the three terms on the right-hand side of equation (1) 
correspond to the stellar ejecta, gas infall, and star formation, 
respectively. ip is the SFR, Yi{m, Z{t — Tm)) is the stellar yield 
of i (in mass units) from a star of mass m and metallicity 
Z{t — Tm), (j^im) is the initial mass function, and Xi is the 
mass fraction of element i. By definition, the sum of Xi over 
all i is unity, and the total surface mass density is identical 
to the integral over the infall rate, niiow and m^p are the 
lower and upper stellar mass limits, respectively, and Tm is 
the main-sequence lifetime of a star of mass m. In practice, 
the first term is split into three equations that deal separately 
with low-mass stars, Type la supernova (SN) progenitors, and 
massive stars. 

We now describe each of the model ingredients in turn, 
distinguishing the inputs that are parameterised using sim- 
ple analytical prescriptions from the nucleosynthesis inputs, 
which are derived numerically from first principles. 

3.1. Analytical Prescriptions 

Initial Mass Function (IMF): The initial mass function de- 
termines the relative birth rate of stars as a function of mass. 
Because different mass stars leave unique chemical signatures 
on the ISM and operate on different characteristic timescales, 
the precise form of the IMF is a key factor driving the evo- 
lution of abundance ratios. This model assumed the Kroupa, 
Tout, & Gilmore (1993) three component IMF. The Kroupa 
et al. function has fewer stars in the low and high mass ends 
of the distribution than the single power-law Salpeter (1955) 
function. In Section 14.2.41 we estimate the sensitivity of the 
results to changes in the upper mass limit of the IMF. 

Star Formation Rate (SFR): A simple analytical law for 
the SFR was adopted. Akin to the Schmidt (1959) law, we 
varied the SFR in proportion to the square of gas surface 
density. The efficiency of star formation was adjusted in order 
to investigate different timescales for metal enrichment. 

Infall Rate: The gas infall rate was assumed to decay expo- 
nentially on timescales ranging between 3 and 9 Gyr. For the 
same final metallicity, a shorter infall timescale corresponds 
to a younger object. The star formation efficiency was also as- 
sumed to be higher in the models that became metal-enriched 
the fastest - this is motivated by the expectation that a deeper 
potential well corresponds to a faster rate of collapse and a 
higher efficiency of star formation. We note that it is the 
combination of both infall rate and star formation efficiency 
that is of foremost importance and these two effects determine 
the time taken to reach a specific metallicity. In the present 
study, the timescale for metal enrichment is more important 
than either the infall timescale or star formation efficiency 
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considered in isolation. For this reason, we present results as 
a function of age and base our discussion on "age-sensitivity" . 

3.2. Stellar Yields and Lifetimes 

Intermediate and low mass stars: We incorporate the 
results of van den Hoek & Groenewegen (1997) whose 
metallicity-dependent yields are specified over the mass range 
0.9-8 Mq. These stars are integral in the synthesis of C and 
N, however observational constraints for these two elements 
are lacking for DLA-B/FJ0812+32. Thus this stellar popula- 
tion will be largely ignored in the following discussions. Low- 
mass stars may be the principle production site for Galactic 
Pb (Travaglio et al. 2001) as well as producing modest quan- 
tities of As, Kr, and Sn via s-process nucleosynthesis. The 
abundance of these s-process elements may prove to be sen- 
sitive indicators of age in high redshift protogalaxies because 
they are generated in 1—4 M© stars with typical lifetimes be- 
tween a few hundred Myr and a few Gyr. However, remaining 
uncertainties in the theory of the s-process, as well as obser- 
vational complications regarding s-element contamination by 
companions stars, may limit the use of Pb as an age indica- 
tor. An investigation into s-process enrichment in DLAs will 
be the focus of a forthcoming paper. 

Type la supernovae (SNe la): A recalculation of the 1986 
W7 model (Thielemann, Nomoto, & Yokoi 1986) by Iwamoto 
et al. (1999) was adopted. About 0.75 M© of iron is ejected 
per event. Our model assumed that 3% of binary systems 
involving intermediate and low mass stars result in SNe la. 
This fraction provides a good fit for the solar neighbourhood 
(e.g. Alibes, Labay, & Canal 2001; Fenner & Gibson 2003). 

Massive stars: Stars more massive than 8 — 10 M0 that 
end their lives in violent supernova explosions are responsible 
for most of the metals in the cosmos. We investigated chem- 
ical evolution using two different sets of SN II yields: 1) the 
Woosley & Weaver grid (1995, hereafter WW95) covering a 
mass range 11-40 M0; and 2) a set of 36 Type II SNe models 
covering the mass range 11-40 Mq at six metallicities (Z = 0, 
10"^, 10"*, 10"^, 6 X 10"^, 2 X 10"^) recently compiled by 
Limongi & Chieffi (2003, hereafter FRANEC) using the latest 
version of the FRANEC code described in Limongi & ChiefB 
(2002). 

For WW95, we took the lower energy "A" models for stars 
< 25 and the higher energy "B" models for heavier stars. 
Taking note of the suggestion by Timmes, Woosley, & Weaver 
(1995) that the WW95 mass cuts may have penetrated too 
deeply within the iron core, we have uniformly halved the iron 
yields from these models ^. The FRANEC models all have 
arbitrary mass cuts corresponding to 0.09-0.1 M© of ^®Ni. 

The nucleosynthesis models only extend to 40 and 35 M0 
for WW95 and FRANEC, respectively. Since stars as heavy 
as 100 Mq are permitted to form in our chemical evolution 
model (albeit ten times less frequently than 40 Mq coun- 
terparts), the predicted yields were linearly extrapolated for 
stars outside the mass grid. Section [4.2.41 discusses the sensi- 
tivity of the results to the assumed upper mass limit on the 
formation of stars. 

Unless otherwise stated, subsequent comparison between 
theoretical and observed elemental abundances will refer to 
the models calculated using WW95. In Section r4.2.5l we show 
that while WW95 and FRANEC make very different pre- 
dictions regarding the magnitude of the "odd-even" effect, 
for most of the elements between O and Ni, for which DLA 
measurements exist, our conclusions drawn from the WW95 
models also apply to the FRANEC models. 

Stellar lifetimes: We adopt metallicity-dependent main- 
sequence (MS) lifetimes calculated by Schaller et al. (1992). 

^ We note that this violates the self-consistency of the stellar 
models, since a shift in mass cut that halves the Fe yield is expected 
to also modify the yield of elements like Co and Ni. 



Stars pollute the interstellar medium with metals over the 
course of their evolution via stellar winds and planetary neb- 
ulae, however this model assumes that all the mass loss takes 
place at the end of the MS phase. This simplification is not 
expected to influence our predictions, since the enrichment 
pattern is mostly controlled by massive stars whose metals 
are returned to the ISM almost entirely during the supernova 
explosion. 

As a word of caution, this type of model assumes that stel- 
lar ejecta is uniformly mixed into the gas phase and therefore 
predicts mean abundance trends. In reality, the products of 
SNe not only mix inhomogeneously with the ambient gas, 
but may also be associated with regions in different ionisa- 
tion states. Despite these caveats, a homogeneous chemical 
evolution model is a reasonable representation of this DLA 
for several reasons: firstly, the measured enrichment pattern 
of this z = 2.626 object reflects the average in the gas phase 
along a line-of-sight; secondly, in order to reach a metallic- 
ity of ~ 1/3 solar, numerous generations of supernovae must 
have exploded, which will tend to "smooth out" the spatial 
distribution of the heavy elements; and, thirdly, similar to 
the majority of DLA (Prochaska 2003) there are only modest 
variations in the relative abundances of the DLA across the 
observed velocity profile. These small variations are likely 
associated to differences in differential depletion implying a 
chemically homogeneous system. 

4. RESULTS 

4.1. Model versus Data 

Figure0presents the difference between the observed DLA 
enrichment pattern and the predicted [X/S]^ from our fiducial 
WW95 model that reaches [0/H] = -0.44 after 2.2 Gyr. The 
dashed line indicates perfect agreement between the model 
and the data, while the dotted lines represent a factor of 
two difference. We conservatively consider predictions falling 
within the dotted lines to be acceptable. The arrows denote 
the direction in which the difference will move if the obser- 
vational limits tighten. For elements like Fe, Ni, and Cr, 
the lower limits are set by dust depletion and are unlikely 
to change. Error bars reflect statistical errors on gas-phase 
abundance measurements, where the statistical error on [X/S] 
equals the quadrature sum of the statistical errors on [X/H] 
and [S/H] given in Tabled When not reported as limits, 
uncertainties in the dust corrections range from from 0.05 to 
0.1 dex, i.e. comparable to the magnitude of the statistical 
errors. 

There is excellent agreement between the fiducial model 
and the DLA data for most elements from B to Ge. This 
agreement is particularly impressive given that the models 
we are applying to a high redshift object were tuned using 
local stellar abundance patterns. Part of the explanation is 
that the DLA pattern is close to solar. Nevertheless, the mod- 
els match an observed trend in the a-elements of decreasing 
[X/S] with increasing atomic number (Figure |5J. Further- 
more, the model reproduces the mildly enhanced "odd-even" 
effect indicated by the limits to the Al, P, and CI abundances. 
There is also good agreement for the iron-peak elements Cr, 
Mn, Fe and Co and the heavier r-process elements Ga and Ge. 
The Ni and Cu predictions do not violate the observational 
limits. The most serious disagreement occurs for Mg and Zn, 

We do not use Fe as the reference element because its nucle- 
osynthesis is subject to many uncertainties relating to the SN II 
mass cut and the evolution of the SN la rate. Oxygen has a more 
robust theory of production, but its observational status and solar 
abundance is less certain. Thus we plot abundances relative to sul- 
fur, for which an ample database of DLA and stellar measurements 
exists. Moreover, the nucleosynthetic origin of S is thought to be 
well understood - being mainly generated by hydrostatic burning 
and explosive O and Si burning in massive stars. 
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Fig. 1. — Logarithmic ratio of the observed and predicted DLA abundance pattern using Woosley & Weaver 1995 SN II yields for a 2.2 Gyr 
model with [O/H] = —0.44. Arrows indicate the direction in which the difference will move if the upper and lower observational limits 
tighten. For elements like Fe, Ni, and Cr, the lower limits are set by dust depletion and are unlikely to change. Error bars reflect statistical 
errors on gas-phase abundance measurements. The level of uncertainty in the dust correction factor is comparable to the statistical errors. 
The dashed line indicates agreement between observations and predictions, while the dotted lines denote deviations of a factor of two. 



which are underpredicted by the model by about 0.45 and 
0.4 dex, respectively. Detailed discussions of these elements 
are found in Sections 14.3.21 and 14.3.51 

4.2. Model Sensitivity 
4.2.1. Dependence on Age 

Figure |5| presents the predicted abundance pattern for a 
model galaxy reaching a final metallicity of [O/H] = —0.44 
after four different periods: 0.46, 0.9, 1.6, and 2.66 Gyr. The 
dust-corrected DLA abundance pattern (as presented in col- 
umn 5 of Table is indicated by black circles, with arrows 
specifying the upper and lower limits, where applicable. Up- 
per limits on As, Kr, Sn, and Pb are shown, although their 
evolution is not modelled. Abundances are on a solar loga- 
rithmic scale relative to sulfur. The dashed line corresponds 
to the solar pattern and the dotted lines are deviations by a 
factor of two. The detection of and limits placed on 23 el- 
ements from B to Pb reveal a roughly solar-like enrichment 
pattern. This may not come as a surprise given the mod- 
erate metallicity of the object. However the redshift of this 
galaxy imposes an upper age limit of ~ 2.5 Gyr, implying a 
much shorter timescale for metal-enrichment than in the solar 
neighbourhood. 

The abundance of most elements relative to sulfur is virtu- 
ally constant over the age range shown in Figure |5| This is 
because the enrichment pattern is mostly set by short-lived 
massive stars. The exceptions are carbon, nitrogen, and the 



iron-peak elements where modest differences are predicted. 
Low and intermediate mass stars have long been held to dom- 
inate the production of N and possibly that of C, although 
theoretical (e.g. Carigi 2000) and observational (e.g. Henry, 
Edmunds, & Koppen 2000) studies also support a strong con- 
tribution to C abundance from Type H SNe. Considerable 
synthesis of C, by the triple-alpha reaction of helium and N, 
via CNO processing, takes place within low and intermedi- 
ate mass stars, which are longer lived than SN II progenitors. 
Hence the enrichment of the interstellar medium with C and 
N is delayed with respect to products of massive star nucle- 
osynthesis. We stress that the results presented in this paper 
are largely insensitive to age. 

The main deviation from the solar pattern is the mild en- 
hancement of the alpha-elements O, Mg and Si with respect 
to Fe or Zn. There is also a trend toward lower relative 
abundances of a-elements with higher atomic numbers. Mas- 
sive, short-lived stars culminating in Type II SNe are the 
primary source of the alpha-elements, while Fe is produced in 
significant amounts by Type la SNe, whose lower mass pro- 
genitors have longer main-sequence lifetimes. Similar alpha- 
enhancement is also seen in halo and thick disk stars where it 
is understood to signify rapid star formation, whereby most 
stars formed from gas that had not yet been polluted with 
iron-rich SN la ejecta. Indeed, the ratio of iron-peak elements 
to sulfur is 30-60% higher in the 2.66 Gyr model compared 
with the 0.46 Gyr model. 
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Fig. 2. — Predicted abundance pattern as a function of age. [O/H] = —0.44 in each model. Abundances are expressed relative to solar 
and scaled to sulfur, such that [X/S] = logio(m{X)/m(S))— logio (m{X)/m(S))Q . Four ages ranging from 0.46 - 2.66 Gyr are represented 
by the symbols indicated. The dust corrected DLA abundance pattern is indicated by lalack circles. Statistical errors on the gas-phase 
abundance are indicated by error bars. Arrows reflect the corresponding lower and upper limits. The dotted lines indicate deviations from 
scaled solar by a factor of two. 



It is apparent from Figure |5| that a protogalaxy that has 
reached 1/3 solar metallicity on a timescale of order 1 Gyr 
has already been contaminated by the nucleosynthetic prod- 
ucts of intermediate mass stars and Type la SNe. This is be- 
cause intermediate mass stars and SN la progenitors can have 
hfetimes as short as a few hundred Myr. Although the rela- 
tive abundances of C, N, and the iron-peak elements exhibit 
moderate age sensitivity in Figure |5| a better probe of age in 
an intermediate redshift protogalaxy might be s-process ele- 
ments such as Kr and Pb, since their origin in low mass stars 
should ensure an enrichment timescale in excess of ~ 1 Gyr 
(Travagho et al. 2001). A preliminary examination of Pb evo- 
lution indicates that [Pb/S] may increase by roughly an order 
of magnitude between 0.5 and 2.5 Gyr (Fenner et al. 2004, in 
preparation). The anticipated detection of s-process elements 
in additional DLAs may prove fruitful in constraining age. 

4.2.2. The Role of Type la SNe 

Figure|5|reveals how the relative contribution from Type la 
SNe increases with the age of the object, for the same metal- 
licity. Type la SNe are believed to be the chief source of 
the iron-peak elements. While the precise identity of SNe la 
progenitors is debated, they must involve intermediate to low 
mass stars in binary systems. The explosion occurs when a 
white dwarf in the binary system accretes enough material 
from its companion to exceed the Chandrasekhar limit. The 
timescale for the SN la event is thus determined by the evo- 



lution of the companion star. The characteristic timescale for 
SNe la is often taken to be 1 Gyr, however this is environment- 
dependent and moreover, the first Type la events may occur 
within only a few hundred Myr (Matteucci & Recchi 2001). 

In theory, the iron-peak abundances might be used to con- 
strain the age of this DLA. This requires assumptions to be 
made about the incidence of SNe la. As mentioned in Sec- 
tion l3.2l we have adopted the Milky Way value for the fraction 
of low and intermediate-mass binary systems that culminate 
in SNe la. In practice, there is a conflict between the predic- 
tions for Mn, which imply an age < 1 Gyr and predictions for 
Fe, which imply that the protogalaxy has been forming stars 
for more than 1 Gyr. In other iron-peak elements, [Co/S] falls 
slightly below the lower limit in all models, while Cr and Ni 
predictions agree with the data. 

From Figure |5| Mn appears to have a very strong depen- 
dence on the age of the simulated galaxy. According to the 
nucleosynthesis prescriptions adopted in this paper, SNe la 
play such a strong role in Mn production that their contri- 
bution exceeds that from SNe II after ~ 2.5 Gyr. It has 
been stated that Timmes et al. (1995) claim Type la SNe 
to be unimportant contributors to Mn synthesis (e.g. Nissen 
et al. 2000), yet it is clear from Figures 4 and 5 from Timmes 
et al. (1995) that Type la SNe produce ~50% of the solar 
Mn abundance. What Timmes et al. showed instead was 
that the same trend of [Mn/Fe] vs [Fe/H] can be obtained 
either with or without SNe la due to the strong metallicity- 
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dependence of the WW95 SN II Mn yields. However, SNe la 
are needed to reach the solar abundance. When the same 
nucleosynthesis prescriptions used in this paper are applied 
to the solar neighbourhood, we find that ~75% of the solar 
Mn abundance originates from SNe la. 

The tendency for Mn to be overpredicted in these DLA 
models may suggest that the main Mn source is metallicity- 
dependent. Other DLA systems have been observed for which 
[Si/Fe] ~ yet [Mn/Fe] is significantly subsolar (e.g. Pet- 
tini et al. 2000). Current nucleosynthesis models cannot 
explain such behaviour, since a solar Si/Fe ratio points to 
a Type la SN contribution, but this is not supported by the 
low Mn abundance, whose origin is understood to lie mostly 
in SNe la. A similar trend has been identified in stellar abun- 
dances of the Sagittarius dSph galaxy (McWilliam, Rich, & 
Smecker-Hane 2003). Invoking metallicity-dependent SNe la 
Mn yields might help explain these apparently incompatible 
observations. 

In a handful of DLAs with low dust content, Pettini et al. 
(2000) discern no correlation between [Mn/Fe] and [Fe/H], 
even though a clear trend is apparent in local stars at all 
metallicities (Carretta et al. 2002). They attribute this to the 
varied histories and morphologies of the objects and suggest 
that a process in addition to Type II and la nucleosynthesis 
may operate to explain Mn enrichment. 

It should be noted that the theoretical yields of Mn from 
WW95 do not have sufficient metallicity-dependence to ac- 
count for the low [Mn/Fe] ratios observed in very metal-poor 
stars (e.g. Alibes, Labay, & Canal 2001). Furthermore, Co/Fe 
ratios up to four times higher than solar are detected in lo- 
cal stars below [Fe/H] ~ —2.5 (Cayrel et al. 2003), in conflict 
with the expectations of SNe II models. Thus, it is likely that 
gaps remain in our understanding of the processes responsible 
for Mn and Co that limit the power of this model to constrain 
the DLA age using these elements. 

4.2.3. Sensitivity to Metallicity 

The production of certain elements (particularly those with 
odd atomic numbers) by massive stars is a strong function of 
the initial stellar chemical composition. We examined the de- 
pendence of the predicted enrichment pattern on the model 
galaxy's metallicity. Figure |21 depicts results from three mod- 
els reaching [O/H] = -0.24,-0.44, and -0.64 after 2.2 Gyr 
of star formation. The ± 0.2 dex variation in [O/H] refiects 
the uncertainty in the measured H I column density. 

The model results displayed in Figure 13 reveal the so- 
called "odd-even" effect, whereby the underabundance of 
odd-numbered elements with respect to their even-numbered 
neighbours increases with decreasing metallicity. This ef- 
fect is a signature of massive SNe from subsolar-metallicity 
stars. Figure |1] shows the production factors relative to S 
from the WW95 SN II models integrated over the Kroupa 
et al. (1993) IMF over the mass range 11-100 Mq. The odd- 
even effect persists from the zero to 1/lOth solar metallicity 
WW95 models, with the trend reversing in the solar model. 
The strength of the odd-even effect in this DLA can be gauged 
by comparing abundances between pairs of neighbouring odd- 
Z and even-Z elements. Evidence for a mild odd-even effect 
comes from the ratios [P/Si] < -0.15, [Mn/Fe] < -0.16, and 
[Ga/Ge] < -0.13. 

Even-elements behave as primary products in SNe II and 
their stellar yields are robust to the initial metallicity. Con- 
versely, the ratios of odd-elements to S show spreads of typi- 
cally ~ 0.2 dex over the 0.4 dex range in [O/H]. The SNe II 
yield of odd-elements from fiuorine to the iron-peak is un- 
derstood to be sensitive to the neutron excess and therefore 
increases with the initial abundance of "seed" nuclei present 
in the star. Thus the ratio of odd-elements to S is highest in 
the [O/H] = -0.24 model. 



Copper exhibits the greatest sensitivity to metallicity, with 
Cu/S having a stronger than linear dependence on O/H over 
the metallicity range considered. Thus [Cu/S] has the poten- 
tial to indirectly constrain the metallicity and the H I column 
density. The upper limit on the Cu detection arises because of 
line blending. However the depletion of Cu onto dust grains 
reduces its abundance in the gas-phase. Accounting for these 
opposing effects, the intrinsic Cu abundance is likely to fall 
within 0.4 dex of the current upper limit. Only the highest 
metallicity models approach the Cu observation. The pre- 
dicted [Zn/S] value, which is underproduced by ~0.4 dex in 
the [O/H] = —0.44 model, is also in closer agreement with 
data for the highest metallicity model. 

Abundances of the odd-elements Al, P, and CI relative to S 
offer a probe into the nature of the "odd-even" effect. How- 
ever, only lower (for Al and CI) and upper limits (for P) exist 
at present in this DLA. These limits are currently not suffi- 
ciently restrictive to eliminate any of the three different metal- 
licity models depicted in Figure 01 Owing to the strength 
of the Al II A 1670 line, only lower limits can be placed on 
the Al abundance in many DLAs, yet absorption-line systems 
present a unique opportunity to acquire empirical evidence for 
the metallicity-dependent behaviour of P and CI, whose weak 
stellar spectral lines mean there is a paucity of data from lo- 
cal stars. A database of measurements of these elements in 
DLAs will help test the theoretical expectation that CI, P, 
and Al follow similar evolutionary paths. 

One might expect from the theoretical predictions shown 
in Figure El that improved observations of P and CI will see 
both [P/S] and [Cl/S] converging toward ~ —0.4. However, 
the abundance of CI shown in these figures was derived from 
the CI I transition, which is not the dominant ion. After cor- 
recting for ionisation, the true CI abundance is likely more 
than 0.5 dex higher than the current lower limit. Even when 
these models are evolved to solar metallicity, Al and CI are 
still underproduced by 0.1 and 0.3 dex, respectively; suggest- 
ing that the odd-even effect may be slightly exaggerated in 
the stellar models, or that intermediate-mass stars provide an 
additional source. It should be noted that the solar reference 
CI abundance used in this paper is ~ 0.15 dex higher than 
the Anders & Grevesse (1989) value. 

Boron is synthesised in a variety of ways including neutrino- 
induced reactions in the shells of SNe II (the r-process) 
and cosmic ray spallation onto seed nuclei in the interstel- 
lar medium (Fields & Olives 1999). The relative importance 
of these processes is expected to vary with metallicity. The 
r-process should produce constant B/O regardless of metal- 
licity, whereas a signature of the cosmic ray process would 
be B/O decreasing with decreasing metallicity. The WW95 
models generate B via the r-process and the satisfactory 
agreement with the observed B/O ratio lends support to a 
r-process origin for B. 

4.2.4. Sensitivity to the Stellar Initial Mass Function 

It has been estimated that uncertainties in the shape and 
limits of the IMF correspond to variations in the absolute 
elemental yields by a factor of 2 (Wang & Silk 1993). We 
have tested the sensitivity of our results to a change in the 
upper mass limit of the IMF from 100 to 40 Mq . Figure ^ 
plots the logarithmic ratio of the abundance pattern obtained 
in these two cases. The abundances of most elements relative 
to S are unchanged. Prominent exceptions are F, Ne, and 
Na, whose yields relative to S decrease by factors ranging 
from 3-5. [O/S] is also lower by ~0.25 dex when the IMF is 
curtailed at 40 Mq . The magnitude of the offset reflects the 
dependence of elemental yield on stellar mass. The elements 
most affected by the shift in upper IMF limit are those whose 
yields increase most steeply with initial stellar mass. The 
possible future detection of F in this DLA may help constrain 
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the IMF. 

4.2.5. Sensitivity to the Type II Supernova 
Nucleosynthesis Prescriptions 

Figure |S| compares the standard WW95 model with the 
same model using FRANEC yields for massive stars. We 
note several main differences: 1) B and F are produced in 
negligible amounts in the FRANEC model because neutrino- 
induced nucleosynthesis processes are omitted from their cal- 
culations; 2) the elements from Cu to Kr are present in trace 
amounts in the FRANEC model; and 3) the odd-even effect 
is much more pronounced in FRANEC models. In particular, 
P, CI, K, and Sc are significantly lower with respect to S in 
the FRANEC model, whereas the abundance of their even-Z 
neighbours matches those from the WW95 model. The deficit 
of CI from the FRANEC models with respect to the obser- 
vations suggests that the odd-even effect may be too severe. 
Despite these differences, there is good consensus between the 
two models for about half of the elements with DLA mea- 
surements and many of the conclusions of this paper can be 
reached irrespective of the choice of SN II yields. 

4.3. Discussion of Individual Elements 
4.3.1. Oxygen 

In the most widely cited set of solar abundances from An- 
ders & Grevesse (1989), the solar oxygen abundance was as- 
sumed to be log(No/NH)0 -f 12 = 8.93. However, recent 
estimates that account for solar granulation and non-LTE ef- 
fects lead to significant downward revisions in the Sun's oxy- 



gen abundance by almost a factor of two (to log(No/NH)0 
-I- 12 = 8.69 - AUende Prieto, Lambert & Asplund 2001). 
Such a shift partly resolves the long-standing dichotomy be- 
tween the Sun's oxygen abundance and that of the local ISM 
(Andre et al. 2003). This work adopts Holweger's (2001) pre- 
ferred value of log(No/NH)0 + 12 = 8.73, whereas WW95 
and FRANEC stellar yields were derived with the Anders & 
Grevesse (1989) values. Owing to uncertainties in stellar nu- 
cleosynthesis inputs such as the ^^C(a,7)^^0 reaction rate, 
convection, and mass-loss (see Langer 1996; Gibson, Loewen- 
stein, & Mushotzky 1997 for useful discussions of these un- 
certainties), we consider our predictions for the O abundance 
to be accurate within a factor of ~ 2. 

4.3.2. Magnesium 

Figure0shows that the [Mg/S] value predicted by the stan- 
dard model is 0.45 dex lower than is observed in this DLA. 
The reason is clear from Figure 2] which shows that Mg/S 
is slightly subsolar in the ejecta from a generation of WW95 
metal-poor stars and tends toward the solar value at higher 
metallicity. The observed value of [Mg/S]DLA = -1-0.3 dex can- 
not be matched by any single generation of WW95 stars and 
nor will any superposition of populations fit the data. Alibes 
et al. (2001) also found that WW95 yields lead to an under- 
production of Mg in their models of the local disk; a prob- 
lem which they suggest could be resolved if intermediate-mass 
stars or SNe la supply additional Mg. Indeed, 4 — 6 Mq stars 
have been shown to contribute to most of the abundance of 
the neutron-rich Mg isotopes at low and intermediate metal- 
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licities (Fenner et al. 2003). 

4.3.3. Sulfur 

Sulfur was used as the reference element in this study in 
preference to Fe because its nucleosynthetic origin is thought 
to be better understood. Hydrostatic burning and explosive 
O and Si burning in massive stars are the main processes 
responsible for S production. Although Fe has traditionally 
been a popular metallicity gauge, owing to easily observable 
stellar spectral lines, the question of its origin is complicated 
by the fact that both Type la and Type II SNe can make 
important contributions to the interstellar Fe content. More- 
over, the Fe yield from massive stars is subject to uncertain- 
ties relating to the SN II mass cut. Oxygen is a common 
alternative reference element, but due to the recent signif- 
icant revisions in the solar O abundance, we chose to plot 
abundance ratios relative to S, which has been detected in a 
reasonable number of DLAs and local stars. 

4.3.4. Titanium 

Titanium is a curious element given that observationally, it 
follows the trends of the alpha-elements in nearby stars, but 
theoretically, is expected to behave as an iron-peak element. 
There is currently only a lower limit on Ti in this protogalaxy, 
however one might expect Ti to be enhanced with respect to 
iron, given the enhancement of O, Mg, and Si. In contrast, 
the model predicts [Ti/Fe] ~ -0.2 dex. The deficit of Ti gen- 
erated by standard SN II and SN la models is a well-known 
unresolved problem (e.g. Timmes et al. 1995, Alibes et al. 
2001). 

4.3.5. Zinc 

Iron is often used as a metallicity gauge in stellar popu- 
lations because of the ease with which it can be observed. 
However, the abundance of iron in gas is difficult to estimate 
because it is readily incorporated into dust grains (Savage & 
Sembach 1996). Furthermore, all of the elements comprising 
the Fe-peak are refractory. Thus Zn has become a popular 
substitute for Fe in DLA studies because it is largely free from 
the effects of dust-depletion and because the Zn/Fe ratio in 
local stars is approximately solar over a wide range of metal- 
licities. The constancy of Zn/Fe vs Fe/H has lead authors to 
propose that Zn, like Fe, is generated mostly by Type la SNe 
(e.g. Matteucci et al. 1993, Mishenina et al. 2002). However 
this need not be the case, since the metallicity dependence of 
Zn production in massive stars can naturally lead to a time 
delay that mimics that of SN la. As can be inferred from 
Figure 1^1 Zn yields predicted by the W7 Type la SNe model 
(Iwamoto 1999) are small in comparison to WW95 SNe II 
yields (whereas the iron-peak elements from V to Ni can be 
largely attributed to SNe la). Figure |3 reveals the strength of 
the metallicity-dependence of Zn production in massive stars. 
These effects can combine to allow Zn and Fe enrichment to 
occur in lockstep despite having different nucleosynthetic ori- 
gins. 

Figure [7| illustrates potential problems associated with us- 
ing Zn as a proxy for Fe. The four panels plot [S/Fe] {left 
panels) and [S/Zn] (rigfit panels) versus time {upper panels) 
and metallicity {lower panels). The curves correspond to dif- 
ferent radii in a Milky Way-like galaxy. The inner radius has 
the strongest star formation rate and the earliest peak. Star 
formation proceeds more slowly with increasing radius. Ow- 
ing to the delayed release of large amounts of Fe from SNe la 
with respect to the release of S from short-lived massive stars, 
the behaviour of [S/Fe] (or [0/Fe]) vs [Fe/H] can be used to 
diagnose star formation histories. Slower and more protracted 
star formation, such as in the outer disk of a galaxy, should 
to lead to subsolar [S/Fe] (and [O/Fe]) at lower [Fe/H] {lower 



left panel). But when Fe is replaced by Zn {lower right panel), 
the behaviour at all radii is virtually indistinguishable. This 
is a consequence of the nucleosynthesis prescriptions adopted 
in this paper that have Zn being produced mostly by massive 
stars as an increasing function of initial stellar metallicity. 
Under this formalism, the evolution of Zn mirrors that of Fe 
only for certain star formation histories, but might deviate 
markedly in other environments. Note that it is impossible to 
attain the solar and subsolar values of S/Zn at the low Zn/H 
that are observed in DLAs (e.g. Centurion et al. 2000) using 
the conventional Type la and II SNe models employed in this 
study, regardless of SF history {lower right panel). This may 
indicate a deficiency in our theories of Zn production from 
the standard SN II and SN la nucleosynthesis models. Other 
authors (e.g. Calura, Matteucci, & Vladilo 2003) have pro- 
duced low [S/Zn] at low metallicity using chemical evolution 
models, but they achieve this by ignoring the theoretically 
predicted SN II and SN la Zn yields and instead assume that 
Zn comes mostly from SNIa and scales with the Fe yield. 

There are two main weaknesses of the standard models of 
Type la and II SNe that highlight the uncertainties in the nu- 
cleosynthetic origins of Zn: firstly, the models are unable to 
explain the supersolar values of [Zn/Fe] in metal-poor stars; 
and secondly, the predicted Zn isotopic composition is at odds 
with solar. Breaking down the theoretical SNe II Zn yield 
into its main isotopes reveals substantial disagreement with 
the solar pattern. The dominant isotope in the Sun is '^^Zn, 
yet WW95 Type II models underestimate its abundance by 
a factor of ~3. Conversely, ^^Zn is overproduced by a similar 
factor. Since the W7 SNe la model also underpredicts ''''Zn, 
other processes might need to be invoked to elevate '^'*Zn to so- 
lar proportions. The uncertainty regarding zinc's production 
sites is unfortunate given its key role in diagnosing DLAs. 

Neutrino-driven neutron star winds (Woosley & Hoffman 
1992; Hoffman et el. 1996) may be an important supplemen- 
tal site for Zn synthesis. WW95 claim that neutrino winds 
accompanying the r-process are probably the chief produc- 
tion site of ^^Zn, as opposed to the neutron capture process 
during helium burning that is included in their SN II mod- 
els. Alternatively, a blend of fallback and mixing in energetic 
hypernovae (HNe) have been found to produce ^*Zn in prodi- 
gious amounts (Umeda & Nomoto 2002), leading to [Zn/Fe] 
values compatible with stellar observations. Each of these 
HNe may eject up to 4x10"'' M0 of ^''Zn; about ten times 
more than a typical SN II. 

Figure [7| illustrates the tendency for a system with a 
"burstier" SFH (such as the innermost disk) to maintain el- 
evated [S/Fe] over a wide range of [Fe/H] because the ISM 
becomes metal-rich before SNe la have a chance to domi- 
nate Fe production. The age-metallicity relation is gradual 
in the outer disk, such that for a given [Fe/H], there has been 
more time for SNe la to contribute to the abundance pat- 
tern. Thus, the outer radii in Figure |7| reaches [S/Fe] = at 
[Fe/H] = -0.7, whereas the inner radii has [S/Fe] = -1-0.2 at 
[Fe/H] = —0.7. Note however, that in the time domain, it is 
the slowly forming systems that maintain elevated [S/Fe] for 
longer temporal periods. Indeed, it takes about 10 Gyr for 
the outer disk in this model to reach solar S/Fe. This corre- 
sponds to a much lower redshift than is characteristic of most 
DLAs. The DLA population is characterised by roughly so- 
lar a/Fe values, a mean metallicity ~ 1/lOth solar, and ages 
< 5 Gyr. According to our standard disk galaxy model, in 
which the gas infall timescale increases with Galactocentric 
radius, only the outer regions of the disk can reach [a/Fe] ~ 
at [Fe/H] ~ —1, but on timescales > 5 Gyr. In order to get 
[a/Fe] ~ at [Fe/H] ~ — 1 on sufficiently short timescales 
one may invoke a short star formation burst of low efficiency 
(to keep the metallicity down) or a galactic wind that prefer- 
entially removes a-elements (e.g. Calura et al. 2003). Dwarf 
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galaxies are more closely associated with these phenomena 
than disk galaxies. 

4.3.6. Gallium and Germanium 

The predictions from our standard model for the abundance 
of Ga and Ge are in rough agreement with the data, however 
we hesitate to draw conclusions based on these results owing 
to the uncertainties surrounding their synthesis. Although 
Woosley & Weaver (1995) include Ga and Ge in their tables 
of yields, they state that the synthesis of all isotopes above 
^^Zn are not considered accurate. Indeed, the WW95 Ga and 
Ge yields are remarkably sensitive to factors including the 
explosion energy. We also note that WW95 predict a stronger 
metallicity-dependence for the even-numbered Ge than for 
the odd- numbered Ga (Figure 2J- 

5. SUMMARY AND FUTURE DIRECTIONS 

We have predicted the chemical evolution of a DLA as a 
function of star formation history, age, and metallicity, and 
find that the enrichment pattern detected in absorption in a 
z = 2.626 protogalaxy is generally consistent with the nucle- 
osynthetic signature of Type II SNe with a moderate contri- 
bution from Type la SNe. Despite this gratifying agreement 
with the data, a few inconsistencies remain that may provide 



insight into nucleosynthetic processes. The underproduction 
of Mg and Zn with respect to S in our models lends sup- 
port to an idea already hinted at by discrepancies between 
the observed abundances in nearby stars and the yields from 
standard SN II models, that these elements may require ad- 
ditional production sites such as: intermediate mass stars, in 
the case of Mg; and neutrino- winds, in the case of Zn. It is 
also likely that tighter constraints on CI and Cu in the future 
will see these elements being underproduced by all but the 
highest metallicity SN II models. This may indicate that the 
odd-even effect is weaker than predicted by stellar nucleosyn- 
thesis models. 

DLA-B/FJ0812-I-32 is expected to soon be complemented 
by a growing database of detailed abundance patterns of 
DLAs covering a range of metallicities and redshifts. These 
observations will represent a new regime for probing chemical 
evolution in diverse environments with an assortment of en- 
richment histories. As well as providing further insight into 
the nature of the DLA population, such a database will com- 
plement local stellar observations to help uniquely constrain 
nucleosynthesis processes. In particular, we await future de- 
tections of CI and P in this and other DLAs. Detection of 
these odd-elements should enable the SN II models of WW95 
and FRANEC to be assessed in terms of their different pre- 
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dictions for the magnitude of the odd-even effect. We also 
propose that the detection of s-process elements like Kr and 
particularly Pb, holds promise as a sensitive measure of the 
age of intermediate to high redshift DLAs. 

Future observations of DLA-B/FJ0812-I-32 are also ex- 
pected to yield measurements of C, N, F, Ga, and Sn. Flu- 
orine, like boron, is difficult to measure in stars, but can be 
detected in absorption in DLAs, providing a test of the r- 
process believed to operate during the core-collapse phase of 
massive stars. It is hoped that Ga measurements will provide 
insight into the r-process in the early universe. Carbon and 
nitrogen detections will be especially fruitful because they can 
be compared against the dataset of N in DLAs and C and N 
in local stars. Furthermore, knowledge of CNO abundances 



in a single DLA will constrain theories of nucleosynthesis in 
intermediate mass stars. 
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